Steering a quantum system over a Schròdinger bridge 
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Abstract 

We outline a new approach to the steering problem 
for quantum systems relying on Nelson's stochastic me- 
chanics and on the theory of Schròdinger bridges. The 
method is illustrated by working out a simple Gaussian 
example. 

1 Introduction 

The problem of steering a quantum system from a given 
state ipo at time ío to another given state V'l at time t\ 
applying an external potential function lies at the heart 
of the field of control for quantum systems. The ap- 
plications are manyfold, including control of molecular 
dynamics, quantum computing, nuclear magnètic reso- 
nance, quantum Communications, etc, see Q-JT^], and 
references therein. In its entirety, the problem consists 
in finding a control potential function {V c (x, t);to < t < 
t±} in a suitable class such that, if Vi(x, t) is the ambient 
potential function, the solution %p{x : t) of the controlled 
Schròdinger equation 



ip{x,t ) = ipo(x), 



V c (x,t)}iP, 



(D 
(2) 



satisfies -0(a;,íi) = ipi(x). Of course, the first natu- 
ral question to pose is the controllability question for 
the Schròdinger equation. Since this is a distributed- 
parameter bilinear system, the problem is highly non- 
trivial, and only very limited results are so far available 
fl. The problem can be effectively studied for N-level 



systems that approximatc in a suitable way the infmitc- 
dimensional system [Q. 



2 Elements of stochastic mechanics and 
Schròdinger bridges 

In this paper, we outline a complctely new approach to 
the steering problem for the Schròdinger equation re- 
lying on Nelson's stochastic mechanic s [ï~^]-pï[ | and on 
the theory of Schròdinger bridges, cf. ]23|, |22[ and refer- 
ences therein. We remark, in passing, that this method 
has nothing to do with the results of [Q-[^6|, where the 
theory of Schròdinger bridges and reciprocal processes 
was employed to construct different versions of stochas- 
tic mechanics. 

Nelson's stochastic mechanics is a quantization proce- 
dure for classical dynamical systems based on diffusion 
processes. Given a quantum evolution {ip(x, t);t < t < 
íi}, where x € M n , a n-dimensional Markov diffusion 
process {x(t);to < t < ti}, called Nelson's process, is 
associated to it as follows. The (forward) Ito differen- 
tial of x(-) is given by 



dx(t) = 



+ 



—V(mogip(x(t),t) + Qlogip(x(t),t)) dt 
m 



dw(t), (3) 



where w is a Standard, n-dimensional Wiener process. 
Moreover, the probability density p(-,t) of x(t) satisfies 



p(x,t) = \^{x,t)\ 2 , Vfc€[ío,*i] 



(4) 



The construction of the Nelson process corresponding to 
ip(x,t) in the case where tj){x,t) vanishes requires con- 
siderable care. It will always be assumed that -0 is of 
class C 2 ' 1 with 



||V^||ieiU[ío,+oo) 



(5) 



This is Carien's finite action condition. Under these 
hypotheses, the Nelson measure may be constructed, 
|Ï5|] , ]Ï7| , Chapter IV], and references therein. Let us 
now outline the theory of the Schròdinger bridges. The 
Kullback-Leibler pseudo-distance between two probabil- 
ity densities p(-) and q(-) is defincd by 



H(p,q) 



log — —p{x)dx. 
q{x) 



This concept can be considerably generalized. Let íl :— 
C([to,ti],M. n ), let W x denote Wiener measure on O start- 
ing at x, and let 

W := I W r dx 



be stationary Wiener measure. Let D be the family 
of distributions on Í7 that are equivalent to W. For 
Q,PgB, we define the relative entropy H(Q, P) of Q 
with respect to P as 



H{Q,P) = E Q [\og 



dQ, 
dP l 



It then follows from Girsanov's theorem that |2£ 
H(Q,P) = H(q(t ),p(t )) + 

rt 



l \[f3l{t)-f} p {t)}-[f}Q{t)-P P {t)]dt 



H{q{t l ),p(t l )) + 
1 



[^{t)-^(t)].[^{t)-^{t))dt 



to 



Here q(to) is the marginal dcnsity of Q at to, [3® and 7^ 
are the forward and the backward drifts of Q, respec- 
tively. Now let po and p\ be two everywhere positivc 
probability densities. Let D(p ,pi) denote the set of 
distributions in D having the prescribed marginal densi- 
ties at to and t\, Given P g D, we consider the following 
problem: 

Minimize H(Q,P) over B(p ,pi). 

This problem is connected through Sanov's theorem 
to a problem of large deviations of the empirical dis- 
tribution, according to Schròdinger original motivation 
p3| , p2[ . If there is at least one Q in D(po,Pi) such 
that H(Q,P) < 00, it may be shown that there exists a 
uniquc minimizer Q* in J}(p 0l pi) called the Schròdinger 
bridge from p to p\ over P . If (the coordinate pro- 
cess under) P is Markovian with forward drift field 



b P .(x,t) and transition density p(a, x, r, y), then Q* is 
also Markovian with forward drift held 

bf(x,t) = bl(x,t) + Vlogcj)(x,t), 

where the function <j) solves together with another func- 
tion 4> the system 



(j)(t, x) = J p(t, x, h,y)<j)(ti,y)dy, 

4>(t, x) = J p(t , y, t, x)(f>(t , y)dy 

with the boundary conditions 

<f>(x, t )cj)(x, t ) = p (x), <j){x, ti)(j)(x, tx) = Pi(x). 
Moreover, we have p(x,t) — (f>(x,t)4>(x,t),Vt e [ío,íi]- 

3 Steering for quantum systems 

We now show that the theory of Schròdinger bridges can 
be employed, jointly with stochastic mechanics, to at- 
tack the steering problem for quantum systems. First of 
all, observe that everything we said about Schròdinger 
bridges continues to hold if we consider finite-energy dif- 
fusions with diffusion coefncient equal to — rather than 
1. Let ipo and V'l be the given initial and final quan- 
tum states. Consider a reference quantum evolution 
{ip{x,t);to < t < ti} solving the Schròdinger equation 

9tp ihj . !„, . , 

and satisfying Carien's finite action condition (||). Let 
P e D be the Markovian measure of the Nelson process 
associated to {ip(x,t)} as in ([|)-(|]). Hence, in particu- 
lar, the probability density satisfies p(x,t) = \4>(x, t)\ 2 . 
Thus, if we write 

i){x, t) = exp[R(x, t) + tS(%, t)], 
the forward drift of the Nelson process is then given by 

òí(x,í) = —VS(x,t) + —VR(x,t). 
m m 

We then have the following result. 

Theorem 1 Let Q* be the Schròdinger bridge from 
\ipo\ 2 to \ipi\ 2 over P (see previous section). Then, Q* 
has forward drift field 

bf(x,t) = — VS(x,t) + —VR(x,t) + — Vlog0(x,í), 
m m m 

where the function <f> solves together with another func- 
tion <j) the system 

^ + (1V5+A V ÍÏ).V^+ A A = O, (6) 
ot m m Zm 



9^ 
~dt 



(— VS- 

m 



-VR)<t> 



2m 



A0 = O,(7) 



with the boundary conditions 

(f>(x,t )$(x,t ) = \tp \ 2 (x), <j)(x,ti)(j){x,ti) = |^i| 2 (x). 
The one-time probability density of Q* satisfi.es 

p(x,t) = (j>(x,t)(f>(x,t). 
Define, for t £ [ío,íi], 

S(M) = S(M) + ft*(M) + £logíp4, (8) 

1 <p{x,t) 

R{x,t) = -log p(x,t). 
Let {ip(x,t);to < t < ti} be defined by 



${x,t) = exp[R(x,t) + -S(x,t)]. 

Then, {^(x^)} solves the Schródinger equation ^) with 
controlling potential function V c (x, t) given by 



V c (x,t) 



V(x)-V t (x) + 

h 2 



&y/p(x,t) _ Ay/p(x,t) 
y/P{x,t) V 'p{x,t) 



(10) 



and we have 

|V>0Mo)| = \^o{x)\, |^(ar,íi)| = \^i(x)\. 

Moreover, the the Schródinger bridge Q* is indeed the 
Nelson process associated to the new quantum evolution 
{$(x,t);t <t <íi}. 

Proof: As is well known [[Ï3|, R and S satisfy the sys- 
tem of nonlinear p.d.e.'s 

noi i 

^ + _viï • VS + —AS = 0, (11) 
dt m 2m 

+ — VS ■ VS* + V - — [Vi? ■ Viï + AR] = 0. (12) 
at 2m 2m 

A long calculation employing definitions (p|)-(p|), and 
equations @, (g), Q and establishes OT. The 
last assertion follows at once observing that (|S|) -([)]) im- 

piy 

bf{x,t) = —VS(x,t) + —VR(x,t). 



Instead, when the ambient potential is zero, we get the 
following remarkable invariance property: 



Fc(M) _^^EI = v{x) ^v^mi 



m y/p(x,t) 



m y/p(x,t) 



The quantity 



h 2 A^p(x,t) 



2m y/p{x,t) 

is called quantum potential in quantum mechanics |Ï4| , 
E7|because it appears in the Hamilton- Jacobi like equa- 



( 9 ) tion (|12|) 



Remark 2 The quantum evolution {ip(x, í);íq < t < 
ti} has, by construction, the desired absolute vàlues at 
times to and t\. Notice that any choice of the "ref- 
erence process" P produces a quantum evolution with 
these properties, and a corresponding control potential 
function. Thus, we can try to choose P (our "free pa- 
rameter"!) so that the phase function S(x,t) has the 
prescribed initial and final vàlues, thereby achieving the 
desired steering, see the example in the following sec- 
tion. 

4 Example 

The aim is that of shifting the mean of a Gaussian wave 
packet, that is, passing from the initial quantum state 



i/j (x) 



U! \ 1 / 4 . U)X 2 , 

-) exp(-— ) 



= exp(Ra{x) + iS (x)) (13) 
at t = to the final quantum state 



Mx) = [-) cxp( ) 

= exp(i?i(x) + iSx(x)) 



(14) 



at t = 1. As we shall show elsewhere pq ], to solve this 
simple problem a direct approach is actually feasible. 
We choose here to solve the problem via the theory of 
Schródinger bridges, in order to illustrate the method 
of the previous section. For notational convenience, let 
us assume that we are in a reference system where m — 
h = 1 and that ui = ir. We will often omit in the sequel 
the function arguments t and x. Let us take as reference 
evolution 



Remark 1 Notice that when V(x) ~ Vi(x), namely the 
reference quantum evolution takes place in the ambient 
potential Vi, the controlling potential reduces to 



V c (x,t) 



A^p{x,t) _ A^p{x,t) 
\/p{x,t) y/p(x,t) 



ip{x, t) = exp(R(x, i) + iS(x, t)) 



where 



S(x,t) = cx + d(t) 



(15) 



R(x,t) = -~( x -m(t)Y (16) 



(17) 



with 



m(t) 



m\ + m 2 í 
m(i) = m 2 



(18) 
(19) 



and d is to be specified. The Schròdinger system (|^)-(0) 
is given by the following two equations for <j> and <f> 



at 



dt 



dt 



dS_ dR\ 90 

dx dx ) dx 

8S_ dR 

dx dx 



l d 2 4> 
20X 1 

2dx^ 



with boundary conditions 



(20) 

(21) 

<t>(x,0)4>(x,0) = p (x) = exp(-wa; 2 ) (22) 

4>(x,l)$(x,Í) = pi(x) =exp{-iü{x-l) 2 ) (23) 

Assuming that 

<j>{x,t) = exp(a(í)x 2 + 0(t)x + j(t)) (24) 

$(x,t) = exp(à(í)a; 2 + 0(t)x +j(t)) , (25) 

we find the following sets of equations for a, à, 0, 0, and 

7,7: 



à - 2aLü + 2a 2 

+ 2ac + 2au}m - cü0 + 2a[í 

b 



j + c0 + oj/3m + \0 2 + a 



and 



à — 2àuj — 2à 2 

+ 2àc + 2àu)m — U10 — 2à0 
7 - lü + c/3 + iü0m - i/3 2 - ce 



(26) 



(27) 



The boundary conditions (p2|)-(p3|) yield the constraints 
on the vàlues of a, à, /3, /3, and 7, 7 at t = 0, 1: 



CtQ + CtQ 

/3o + /3o 
7o + 7o 



= — Lü 

= 
= 



and 



Q'l + a\ = 



(28) 



(29) 



(30) 
(31) 
(32) 



satisfy The value of 0{t) can be found by 

integration 

^(í) = 2wm(í)-e- art (/3 +2wm 1 ). (33) 



/3i +/3i = 2w 

7l + 7l = ~U) 

It is easy to see that 

a(t) = 
à(t) = —Lü 
0{t) = /3 e wt 



Using the relation 0i + $1 — 2u>, we obtain 

1 + m>ie~ u - (m\ + m 2 ) 
00 = 2üj 



(34) 



Given /3 and /3, it is possible to obtain 7 and 7 (by 
integration), where 70 is so far to be specified. Thus, we 
get 



7(*) 
7(í) 



4ui ' 



70 + HL(1 - e 2 - 4 ) + AK ~ e" 4 m(í)l35) 



-7o + ^[/3 2 (l-e- 2wt ) 

+4/3 we- 2wt (e tüt m(í) - 
+4iü 2 e- 2uJt {e UJt m(t) - mi) 2 ] 



(36) 



It can be checked that the constraint 71+71 = — w is 
always satisficd, whatever the value of 70. Matching of 
the phase at t — 0, 1 requires that 

S(x,0) +R(x,0) +ln <j){x,0) = S Q (x) + -ln p {x) 

S(x,l) +R(x, 1) +ln <f>(x,l) = Si(x) + -ln p x {x) , 

or equivalently that 

m 2 + ljtoi + 0o = (37) 
7o _| m 2 + do = (38) 

m 2 + w(mi + m 2 ) + /3i - lü = (39) 
|(mi+m 2 ) 2 +7 1 + |+c?i = 0, (40) 

where do and <ii are the vàlues of d{t) at í = 0, 1. Some 
computations show that eqns. ( ^4|) and (|37j)-(40) are 
satisfied by the following vàlues 



1 



m 2 

0o 

7o 
do 



2 - 2e u + uj + we w 
u(e" + 1) 
2 - 2e" + w + oje w 
2ui 

2 - 2e u + lü + we" 
u) (e" - l) 2 

2 (2 - 2e"> + üj + we") 2 



= 



di = - 



Lü(l + e w )(l - e" + uj + ue u ) 



(2 - 2e u + lü + iüe u 
A possible choice for d(t) is then 
d(t) = dit . 



(41) 

(42) 

(43) 

(44) 
(45) 
(46) 

(47) 



By choosing m, d, a, à, 0, 0, and 7, 7 according to equa- 
tions both the reference 
evolutio n (^5[ ) and the solution ((/>, (j)) of the Schròdinger 
system ([20|)-(p3"|) are completely specified. The expres- 
sions of the controlled quantum evolution ip(x,t) and of 
the controlling potential function V c {x,t) can then be 
derived as in Theorem 1. 
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